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We investigate the quark mass dependence of baryon masses in 2+1 flavor lattice QCD using 
SU(3) heavy baryon chiral perturbation theory up to one- loop order. The baryon mass data used 
for the analyses are obtained for the degenerate up-down quark mass of 3 MeV to 24 MeV and 
two choices of the strange quark mass around the physical value. We find that the SU(3) chiral 
expansion fails to describe both the octet and the decuplet baryon data if phenomenological values 
are employed for the meson-baryon couplings. The SU(2) case is also examined for the nucleon. 
We observe that higher order terms are controlled only around the physical point. We also evaluate 
finite size effects using SU(3) heavy baryon chiral perturbation theory, finding small values of order 
1% even at the physical point. 

PACS numbers: 11.15.Ha, 12.38.-t, 12.38.Gc 



I. INTRODUCTION 

The aim of the PACS-CS project is full QCD calcula- 
tions on the physical point avoiding any contamination 
due to chiral extrapolations. At the first stage of this 
project [l| we have succeeded in reducing the up-down 
quark mass from 67 MeV, the minimum value reached by 
the previous CP-PACS/JLQCD work@, to 3 MeV, cor- 
responding to the decrease of pion mass from 702 MeV to 
156 MeV. This work allowed us to make detailed chiral 
analyses on the pseudoscalar meson sector with the use of 
chiral perturbation theory (ChPT). An important find- 
ing is that the strange quark mass is not small enough to 
be treated by the SU(3) ChPT up to the next-to-leading 
order (NLO). For the octet and decuplet baryon masses 
we simply employed a linear chiral extrapolation to the 
physical point assuming isospin symmetry and analytic- 
ity of the strange quark contribution around its physical 



SU(2) BChPT analyses in 2 flavor lattice QCD were pre- 
viously made by other collaborations and reported in 
Refs. [H, We also discuss the magnitude of finite 

size effects based on the NLO SU(3) HBChPT. 

This paper is organized as follows. In Sec. [IT] we briefly 
review the results obtained in Ref. [l[ to make the pa- 
per self-contained. In Sec. HID we apply SU(3) HBChPT 
analyses to the octet baryon masses. We present the fit 
results for the LECs and discuss the convergence behav- 
ior up to NLO. The same analysis is repeated for the 
decuplet baryon masses. Section [TVl describes the fit re- 
sults of the nucleon mass with the SU(2) HBChPT. In 
Sec. [V]we discuss the magnitude of finite size effects for 
the baryon masses based on the SU(3) HBChPT. Our 
conclusions are summarized in Sec. I VI I 



II. SU(2) LINEAR CHIRAL FIT FOR BARYON 
MASSES 



In this article we reinvestigate the quark mass depen- 
dence of the octet and decuplet baryon masses employ- 
ing the SU(3) heavy baryon chiral perturbation theory 
(HBChPT) up to NLO 0, Hi]. The results are compared 
with those of the linear chiral extrapolation obtained in 
Ref. [l[. We also examine the convergence property of 
HBChPT fits to the lattice results. For the nucleon mass 
we examine the SU(2) HBChPT with an analytic ex- 
pansion of the low energy constants (LECs) in terms of 
the strange quark mass around its physical value. The 
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We give a quick review of the chiral analyses employed 
in Ref. [l|. The extrapolation to the physical point is 
performed with the following fit formula: 



/3m 



AWI 
ud 



■ 7 TO S 



AWI 



(1) 



where denotes the Axial Ward Identity quark mass 

for the up-down quark and m^ WI for the strange quark. 
Since we choose k s around the physical strange quark 
mass, the above formula is essentially an SU(2) chiral 
expansion with the strange quark contribution analyti- 
cally expanded around its physical value. Actually we 
can rewrite the formula UJ as 
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where a' — a + jmf^ with the physical value of 

the strange quark mass. This is a linear expansion of the 
baryon mass around (ra^ 1 , m^ WI ) = (0, m^Sf). 

The simulation points and the measured hadron 
masses are given in Table I and Table III in Ref. [l|, 
respectively. The fit range is the lightest four points at 
«ud > 0.13754, where the pion mass varies from 156 MeV 
to 410 MeV and m^f 1 from 3 MeV to 24 MeV in the MS 
scheme. In Figs. Q] and [2] we present the fit results for the 
octet and the decuplet baryon masses, respectively. Star 
symbols denote the extrapolated values at the physical 
point whose numerical values are listed in Tables HI and HT1 
The physical point together with the lattice cutoff is de- 
termined with ttitt , rriK , tuq inputs by applying the NLO 
SU(2) ChPT fit to the PS meson sectorpj. The data are 
reasonably described by the linear function Eq. (TT]) . The 
values for a, j3, 7 and % 2 /dof are summarized in Tables 
IIIII and HVl In oder to investigate the convergence behav- 
ior the contribution of each term in Eq. Q is drawn in 
Fig. [3] for the octet baryon masses, where m^ WI is fixed 
at the measured value for (k u( j, k s ) = (0.13754,0.13640). 



We observe that the 0(m s — m s p h) terms are small, and 
that 0(m u( j) contributions are less than 20% of that of 
a 1 for m^ l <0.01. This suggests that higher order terms 
in m u( j and m s — m SjP h would be small. Similar trends 
are found for the decuplet baryon masses in Fig. |4j 



III. SU(3) HBCHPT ANALYSES ON BARYON 
MASSES 

A. Lagrangian to leading order 

We use the continuum HBChPT in this article, leaving 
the development of the Wilson HBChPT, which incorpo- 
rates the chiral symmetry breaking effects of the Wilson- 
type quark action, for future work. We follow the nota- 
tion described in Ref. [5j . The lagrangian of HBChPT is 
written in terms of velocity-dependent baryon fields with 
a perturbative derivative expansion. To fix the notation 
we present the leading order terms: 



C = (Biv ■ DB) + 2a M (BBM+) + 2f3 M (BM+B) + 2a M (BB) tr (M. 
- (f» {iv-D- A} T M ) + 2 7M (T^M+T^) - 2a M tr (M+) 

+2a (BS^BA,,) + 2/3 (BS^A^B) + 2H (T , S»A lt T u ) 
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where B and T represent the velocity-dependent octet 
and decuplet baryon fields with the four-velocity 
u M , respectively. This lagrangian contains 9 LECs: 
O-m, Pm 1 CM) 7m j a i Pi T~L) C- Hereafter we use the ax- 
ial couplings F, D instead of a, /3. They are related with 



2F, 
+ F. 



(4) 
(5) 



The decuplet-octet mass difference denoted by A = 
hit — m-B has a comparable magnitude to inter- multiplet 
mass splittings of both the octet and the decuplet. We 
are allowed to treat it as a perturbation. The octet pseu- 
doscalar mesons couple derivatively to the baryon fields 
through the vector combination V^, the axial vector one 
A^ and the chiral covariant derivative D^. The light 
quark masses m u , rn&, m s are contained in We refer 

to Ref. ,5| for the explicit expression of V^, Ap, D^,M + . 

In this section the physical to u( j and m s are determined 
by the SU(3) ChPT analyses in the pseudoscalar meson 
sector including and rriK', tuq is an additional input 
for the lattice cutoff 111. 



B. Octet baryons 

In the SU(3) HBChPT the renormalized mass of the 
z-th octet baryon is expanded as 



(1) (3/2) 

rriBi = rn-B — m B ~ m B 



(6) 



where is the octet baryon mass in the chiral limit un- 



der the SU(3) flavor symmetry, and is the 0(m^ t; ) 
contribution to the i-th octet baryon. The LO correc- 
tions are given by 
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with m u( j the averaged up-down quark mass, m s the 
strange quark mass, and ajf, Pm, &m are LECs in the 
lagrangian (|3|). 



The NLO contributions, which are obtained from one- 
loop diagrams, are written as 
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where A is the octet-decuplet baryon mass difference in defined in the SU(3) chiral limit with the convention of 
the SU(3) chiral limit. The pion decay constant / is also /„. = 130.4 Me V. The function T is expressed as 
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with T(m,p, 0, /it) = 7rm5, and /i the rcnormalization 
scale. This formula assumes A > m^. For A < 
we apply the analytic continuation: 



A 2 — m\ + iehi 



A 2 — m 2 + ie 
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We revisit the function T later in Sec. [V] to discuss fi- 
nite size effects. The contributions from the octet-octet- 
and the decuplet-octet-axial couplings are factored out 
by A?* and C 2 B?\ respectively. We summarize their 
values in Table [Vj The LECs are phenomenologically 
estimated as@, \ 



D = 0.80, F = 0.47, C=1.5. 



(14) 



We first present the fit results up to leading order 
employing the formula = tiib — m^ . The values 
for ms, ajw, /3m, om are given in Table IVll together with 
% 2 /dof. We also present the extrapolated values at the 
physical point in Table[T]in comparison with the results of 
the SU(2) linear fit Eq. ([T]). The two set of results are con- 
sistent within 2a error. Figure [5] shows that the data are 
reasonably described by the formula. The convergence 
behavior in Fig. [SJ however, is disappointing because of 
the sizable contribution of the 0(m s ) term. This point 
is in a sharp contrast to the SU(2) linear expansion of 
SecHJ where the contributions of the 0(m u d, m s — m s>p h) 
terms are well controlled in the range of to u( j<0.01. 

In the NLO fit the number of LECs increases up to 
seven. We give the fit results in Fig. [7] and Tables IVT1 and 
UJ Although the fit works in a reasonable manner, a crit- 



ical observation is that the results for D,F,C are essen- 
tially consistent with zero showing a significant deviation 
from the phenomenological estimates in Eq. (fT4|) . To ex- 
amine the contributions of the NLO terms, we make a 
fit with £), F, C fixed at the phenomenological estimates. 
Figure [8] and Table IVII show that this fit assumption is 
strongly disfavored because of a prohibitively large value 
of x 2 /dof. The reason is found in Fig. O If D,F,C are 
fixed at the phenomenological estimates, the magnitude 
of the NLO contribution is 2—5 times larger than that 
of the data even at m u( j = 0. We also remark that the 
same situation holds if either of D,F,C is fixed at the 
phenomenological estimate. 



C. Decuplet baryons 

Let us turn to the decuplet baryons. The chiral expan- 
sion of the i-th decuplet baryon mass is written as 



rriTi = nriT + vrirp 



(i) 



,(3/2) 



(15) 



with ttit the decuplet baryon mass in the SU(3) chiral 
limit. The LO and NLO corrections are given by 
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]T -H 2 A^(m , 0, M ) + C 2 Bj^(m , -A, M ) 



(20) 



where 7m,cm,H are additional LECs, and TL 2 A^ 

and C 2 B^ denote the contributions coming from the 
decuplet-decuplct- and the dccuplet-octet-axial cou- 
plings, respectively, whose complete list is given in Ta- 
ble [Vn] Note that the function T(m,^,— A,/i) yields an 
imaginary part if A > m0[5(. Our analyses are made 
by taking account of the real part only. In the region 
of A > ?Ti0 it might be problematic to apply the contin- 
uum chiral expansion for infinite spatial volume to lattice 
data. We leave the proper treatment to future work. 

The LO fit results with the formula = tut + to^P 
are presented in Figs. ITTJI and [TT1 and Tables [Villi and [Hi 
The situation is similar to the case of the octet baryon 
masses. The LO formula successfully describes the quark 
mass dependence of the octet baryon masses. The ex- 
trapolated values at the physical point do not show any 
sizable deviations beyond error bars from those of the 
SU(2) linear fit Eq. |T]). The contributions of the 0(m s ) 
term, however, are rather large compared to the magni- 
tude of tut- The SU(2) linear formula shows a better 
convergence behavior than the SU(3)-symmetric linear 
formula. 

The NLO fit is carried out with and without fixing C at 
the phenomenological estimate. The results for the lat- 
ter case are given in Fig. [12] and Tables [Villi and HT1 We 
find a reasonable value for x 2 /dof, though the vanishing 
result for C shows inconsistency with the phenomenologi- 
cal estimate. Figure[T3l illustrates the situation whenC is 
fixed to the phenomenological value: The NLO fit hardly 
describes the quark mass dependence, and suffers from 
an unacceptable value of x 2 /dof. In Fig.Q3]we find that 
the convergence is worsened by the NLO contributions. 



IV. SU(2) HBCHPT ANALYSES ON NUCLEON 
MASS 



of m s — TO s ,ph yields higher order corrections beyond 

We apply the formula (|2"Tj) to the four data points with 
O77i u d^0.01, employing the experimental axial coupling 
cja = 1.267, and the value of / already determined from 
the SU(2) ChPT fit for m„, f„,f K in Ref. Q. The re- 
sults are given in Fig. [15] and Table IIXI The value of 
X 2 /dof is sufficiently small. However, we should remark 
two points. Firstly, Fig. [15] shows that the fit results fail 
to predict the quark mass dependence of the data be- 
yond to^ 1 = 0.01. The extrapolated value at the phys- 
ical point, for which we obtain m^v = 0.382(25), also 
undershoots the experimental one sizably. Secondly, the 
LO and NLO contributions quickly increase as the quark 
mass increases so that the good convergence region is 
restricted near the chiral limit. 

It is worthwhile to make a comparison of our fit 
results with the 2 flavor twisted mass case given in 
Ref. [l(| where the authors apply the NLO formula 
(f2Tj) to the nucleon mass choosing gA = 1.2695(29) and 
/„ = 0.092419(7) (25) GeV[l8j. A particular interest 
exists in the parameter c\ which is responsible for the 
pion-nucleon sigma term: a(t = 0) = — 4cito 2 with t 
the squared momentum transfer to the nucleon. Our re- 
sult is ci = -1.44(31) GeV" 1 with a" 1 = 2.176 GeV, 
which is consistent with ci = —1.19(1) GeV in the con- 
tinuum limit of the 2 flavor twisted mass case. We ob- 
tain cr(0) = 75(15) MeV for the sigma term, which prefers 
cr(0) = 64±7 MeV given by a recent analysis with the new 
experimental data[ll| to the conventional phenomenolog- 
ical estimate cr(0) ~ 45 MeVpjfl. A reliable way for the 
precise determination of the sigma term, however, is the 
direct calculation of the forward matrix element of the 
nucleon (N\qq\ N) which requires both the connected and 
the disconnected diagrams 



quires both 

51 El El. 



In the framework of SU(2) HBChPT the nucleon mass 



up to 0(m^J ) is given by 



A ? ^9 A s 

m N = m Q - Ac x m n - „ bi;, 

167T/^ 



(21) 



where too and / are defined in the SU(2) chiral limit, C\ 
is a low energy constant and gA denotes the axial vector 
coupling of the nucleon. We further expand mo around 
the physical strange quark mass: 



m = m + m' (m s - m SiPh ). 



(22) 



Since \m s — 777, SiP h| is comparable to to u( j in our sim- 
ulations, the analytic expansion of c\,gA,f in terms 



V. FINITE SIZE EFFECTS 

The one-loop correction for the baryon mass is evalu- 
ated by the following function, 



-8tt' 



d 4 p 



P 2 



o (27r)Mw-A)(^+^ + m*)' 



(23) 



where the integral is defined in the Euclidean space-time. 
This leads to Eq. (fT2"|) after the dimensional regulariza- 
tion and the renormalization in the MS scheme with the 
scale [i. In a finite spacial volume of linear size L, the 
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integral over the spatial components of the loop mo- 
mentum p is replaced by a sum over discrete momenta 
p = (2n/L)n, 



where we assume that the time direction is infinite. We 
define the finite size correction for T as 



•F (L) (m ,A) 
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where /3% = A 2 + 2AA + mi and 
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With the use of the master formula! 16, 11 
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where K n (z) is a modified Bessel function of the sec- baryon masses at next-to-leading order are 
ond kind, one finds that the finite size corrections to the 
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In order to evaluate the magnitude of the finite size 
effects, let us consider the asymptotic form of 5l3~ 17]: 



= — 6\/27rf7j 
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The finite size effect for the baryon mass <5lti = m(L) — 
m(oo) is expressed as 



= A 
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where we neglect the sub-leading contributions from tuk 
and m^. An intriguing point is that E\ and E 2 diminish 
as the pion mass decreases if the product of m v L is kept 
fixed. The values of E\ and E 2 with ah — 32 at the 
physical point can be evaluated as 



aE 1 
aE 2 



6.61 x 1CT 



1.43 x 1CT 4 , 

where we employ the following results in Ref . [l| : 

m ld B a = 0.00859(11) [GeV 2 ], 
af = 0.0546(39), 
a' 1 = 2.176(31) [GeV] . 



(36) 
(37) 



(38) 
(39) 
(40) 



We also use A = mr — mu — 0.16 obtained by the LO 
SU(3) HBChPT fit for the baryon masses. Table H sum- 
marizes the coefficients A, B and the relative finite size 
correction normalized by the baryon mass extrapolated 
at the physical point with the formula ((T|). To evaluate 
the numerical value of A and B we use the phenomeno- 
logical estimates D = 0.80, F = 0.47 and C = 1.5 and 
assume Ti = 1 which is comparable to our fit results in 
Table lyTITl We find that the HBChPT predicts fairly 
small finite size corrections even at the physical point: 



The magnitude is less than 1% for all the channels. We 
additionally remark that the f2 baryon mass is free from 
the leading contribution of the finite size effects. This is 
another fascinating feature to choose mo as one of the 
physical inputs. 



VI. CONCLUSION 

We have investigated the chiral behavior of the octet 
and the decuplet baryon masses based on the SU(3) 
HBChPT. At LO we find reasonable fit results both for 
the octet and the decuplet baryon masses, though rather 
large strange quark contributions are observed. This 
point is contrary to the SU(2) linear chiral expansion 
where the LO strange quark contribution is well con- 
trolled around the physical m s . Inclusion of the NLO 
contributions makes the situation worse: The fit results 
are incompatible with the phenomenological estimates 
for the low energy constants D,F,C both for the octet 
and the decuplet cases. We have also applied the NLO 
SU(2) HBChPT to the nucleon mass. The quark mass de- 
pendence is reasonably described below m^^O.Ol. The 
good convergence property, however, is observed only 
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near the chiral limit. The finite size effects predicted 
by the SU(3) HBChPT turn out to be fairly small: at 
most less than 1%. However, we should be aware of the 
possibility that the value m n L w 2 at the physical point 
for the aL — 32 lattice may be beyond the applicabil- 
ity range of Eq. (f35|) . Comparisons with lattice data on 
larger lattices should settle the issue here. 

In order to avoid the difficulties associated with the 
chiral extrapolation of the baryonic quantities we are now 
carrying out a simulation directly on the physical point. 
This simulation is being made on a larger lattice size so 
that a direct study of finite size effects is possible. 
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TABLE I: Octet baryon masses at the physical point. NLO results are obtained without (case 1) and with (case 2) fixing 
D, F,C at the phenomenological estimate. 





linear 


LO 


case 1 


NLO 


case 2 


N 


0.438(20) 


0.4532(75) 


0.447(15) 




0.322(32) 


A 


0.502(10) 


0.5199(72) 


0.517(11) 




0.387(82) 


E 


0.531(11) 


0.5439(77) 


0.546(12) 




0.53(16) 




0.5992(75) 


0.5987(80) 


0.600(14) 




0.62(18) 



TABLE II: Decuplet baryon masses at the physical point. NLO results are obtained without (case 1) and with (case 2) fixing 
C at the phenomenological estimate. 





linear 


LO 


case 1 


NLO 


case 2 


A 


0.586(19) 


0.612(9) 


0.604(19) 




0.545(20) 


E* 


0.657(15) 


0.666(9) 


0.663(14) 




0.604(17) 




0.718(12) 


0.720(9) 


0.721(14) 




0.694(17) 


Q 


0.769(11) 


0.774(11) 


0.777(17) 




0.815(21) 



TABLE III: Fit results with the linear formula Eq. ([T]) for the octet baryon masses. 





N 


A 


E 




a 


0.371(50) 


0.375(28) 


0.381(33) 


0.408(21) 


P 


11.6(2.4) 


9.6(1.0) 


8.0(1.0) 


4.24(46) 


7 


1.8(1.3) 


3.90(76) 


4.7(9) 


6.24(58) 



X 2 / dof 0.63(2.5) 1.3(2.3) 0.6(1.4) 0.09(59) 



TABLE IV: Same as Table [TTTI for the decuplet baryon masses. 





A 


E* 


— '4-' 


Q. 


Q 


0.527(60) 


0.538(48) 


0.548(40) 


0.552(33) 


P 


10.7(2.1) 


6.3(1.5) 


3.41(82) 


1.80(57) 


7 


1.6(1.6) 


3.8(1.3) 


5.5(1.1) 


7.2(9) 



X 2 / dof 0.4(1.4) 0.2(1.3) 0.00(21) 0.6(1.7) 



TABLE V: Coefficients for the octet baryons A~? 1 and B^' . 













<t> 






7T 


K 


V 


7T 


K 


V 


N 


U D + F f 


|(5D 2 -6DF + 9F 2 ) 


ICD-3F) 2 


4 
3 


i 

3 





A 


2D 2 


l(D 2 +9F 2 ) 


ID 2 


1 


2 
3 





E 


l(D 2 + 6F 2 ) 


2(D 2 + F 2 ) 


§£> 2 


2 
9 


10 
9 


1 

3 




l(D~F) 2 


|(5D 2 + GDF + F 2 ) 


l(D + 3F) 2 


1 

3 


1 


1 

3 
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TABLE VI: Fit results with the SU(3) HBChPT for the octet baryon masses. NLO results are obtained with (case 2) and 
without (case 1) fixing D,F,C at the phenomenological estimate. 





LO 


NLO 

case 1 


case 2 


phenom. 


tub 


0.410(14) 


0.391(39) 


-0.15(9) 




OLM 


-2.262(62) 


-2.62(62) 


-15.3(2.0) 




Pm 


-1.740(58) 


-2.6(1.5) 


-21.3(3.0) 




om 


-0.53(12) 


-0.71(34) 


-9.6(1.4) 




D 




0.000(16) xl0~ 8 


0.80 fixed 


0.80 


F 




0.000(9) xl0~ 8 


0.47 fixed 


0.47 


C 




0.36(30) 


1.5 fixed 


1.5 



X 2 /dof 1.10(63) 1.39(77) 153(82) 



TABLE VII: Coefficients for the decuplet baryons AT* and . 







A Ti 
A 4> 












TV 
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5 
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E* 


4 
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!) 
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i — i* 
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1 
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6 


3 


3 


3 








2 
3 


2 
3 





4 
3 






TABLE VIII: Fit results with the SU(3) HBChPT for the decuplet baryon masses. NLO results are obtained with (case 2) and 
without (case 1) fixing D, F, C at the phenomenological estimate. 





LO 


case 1 


NLO 


case 2 


phenom. 




0.570(16) 


0.550(43) 




0.359(43) 




JM 


2.745(80) 


3.4(1.1) 




3.88(25) 




0~M 


-0.56(15) 


-0.96(75) 




-2.88(40) 




c 




0.000(21) xl0~ 8 




1.5 fixed 


1.5 


n 




0.50(49) 




0.000(4) xlO -8 





\ ' dof 0.46(48) 0.50(60) ~ 21.5(9.5) 



TABLE IX: Fit results with SU(2) HBChPT up to NLO for the nucleon mass. 



NLO 


mo 


0.258(63) 


m' 


2.8(1.4) 


Cl 


-3.14(68) 


X 2 /dof 


0.2(9) 
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TABLE X: Relative finite size correction normalized by the baryon mass extrapolated at the physical point with the formula 
Eq. JT]). LECs are chosen to be D = 0.80, F = 0.47, C = 1.5 and H = 1. 



A 



B 



3(D + F) 2 
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2C 2 
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c 2 



J? = fcm/m[%] 
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0.45 
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FIG. 1: Fit results with the linear formula Eq. ([TJ for the octet baryon masses. Experimental values are given in lattice units 
with a" 1 = 2.176 GeV in Ref. [l[. 
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FIG. 2: Fit results with the linear formula Eq. |T} for the decuplet baryon masses. Experimental values are given in lattice 
units with cT 1 = 2.176 GeV in Ref. Q. 



14 




0.01 0.02 
AWI 

m ud 



0.03 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0. 







m . 



0.01 0.02 
AWI 

m u d 



0.03 




AWI 

m ud 



0.03 



0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0, 







0.01 



0.02 



AWI 

m u d 



0.03 



FIG. 3: 

at (k uc i 
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FIG. 5: Fit results with the SU(3) HBChPT up to LO for the octet baryon masses. Experimental values are given in lattice 
units with a" 1 = 2.176 GeV in Ref. [1. 
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FIG. 6: Convergence behavior for the octet masses with the SU(3) HBChPT up to LO. m s is fixed with the measured value at 
(«ud, «s) = (0.13754, 0.13640). 
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FIG. 7: Fit results with the SU(3) HBChPT up to NLO for the octet baryon masses. Experimental values are given in lattice 
units with a" 1 = 2.176 GeV in Ref. [1. 
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FIG. 8: Fit results with the SU(3) HBChPT up to NLO for the octet baryon masses. D, F,C are fixed at the phenomenological 
estimates. Experimental values are given in lattice units with a -1 = 2.176 GeV in Ref. [lj]. 
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FIG. 10: Fit results with the SU(3) HBChPT up to LO for the decuplet baryon masses. Experimental values are given in 
lattice units with a' 1 = 2.176 GeV in Ref. Q. 
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FIG. 11: Convergence behavior for the decuplet masses with the SU(3) HBChPT up to LO. m s is fixed with the measured 
value at (K ud ,K s ) = (0.13754,0.13640). 
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FIG. 12: Fit results with the SU(3) HBChPT up to NLO for the decuplet baryon masses. Experimental values are given in 
lattice units with a -1 = 2.176 GeV in Ref. [J- 
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FIG. 13: Fit results with the SU(3) HBChPT up to NLO for the decuplet baryon masses. C is fixed at the phenomenological 
estimates. Experimental values are given in lattice units with a" 1 — 2.176 GeV in Ref. 0|. 
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FIG. 14: Convergence behavior for the decuplet masses with the SU(3) HBChPT up to NLO. C is fixed at the phenomenological 
estimates. m s is fixed with the measured value at (ft u d, k s ) = (0.13754, 0.13640). 
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results with the SU(2) HBChPT up to NLO for the nucleon mass (left) and comparison of LO and NLO 
(right). Experimental values are given in lattice units with a -1 = 2.176 GeV in Ref. 



